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Motivated by experiments on NiGa2S4, we discuss characteristic (finite temperature) properties of spin S = 1 
quantum antiferromagnets on the triangular lattice. Several recent theoretical studies have suggested the possi- 
bility of quadrupolar (spin-nematic) ground states in the presence of sufficient biquadratic exchange. We argue 
that quadrupolar correlations are substantially more robust than the spin-nematic ground state, and give rise to 
a two peak structure of the specific heat. We characterize this behavior by a novel T > semiclassical approx- 
imation, which is amenable to efficient Monte Carlo simulations. Turning to low temperatures, we consider the 
effects of weak disorder on incommensurate magnetic order, which is present when interactions beyond nearest 
neighbor exchange are substantial. We show that non-magnetic impurities act as random fields on a component 
of the order parameter, leading to the disruption of long-range magnetic order even when the defects are arbi- 
trarily weak. Instead, a gradual freezing phenomena is expected on lowering the temperature, with no sharp 
transition but a rapid slowing of dynamics and the development of substantial spin-glass-like correlations. We 
discuss these observations in relation to measurements of NiGa2S4. 

PACS numbers: 75.10.Jm,75.10.Hk 



I. INTRODUCTION 

The triangular lattice antiferromagnet has long been of in- 
terest because of its potential to exhibit exotic phases as a re- 
sult of frustration and because it is the underlying lattice of 
many real materials. Unfortunately though, the simplest tri- 
angular lattice models have a well known tendency to order 
into states that are basically classical. For instance, it is gen- 
erally believed that the nearest-neighbor spin 1/2 triangular 
lattice antiferromagnet orders into the 120-degree spiral state 
at zero temperature. However, since the models that describe 
real triangular lattice systems generally contain additional in- 
teractions, the possibilities for more interesting behavior re- 
main vast. 

One triangular lattice system that was recently proposed as 
a candidate for exotic behavior was the spin one antiferromag- 
net NiGa2S4, which exhibits no long-range spin ordering at 
zero temperature despite a low temperature susceptibility and 
specific heat' that suggest two-dimensional magnons; in par- 
ticular, the specific heat grows very clearly as below lOA'. 
The specific heat furthermore shows no signs of a phase tran- 
sition but instead shows two rounded peaks at lOK and 80K, 
the latter being the Curie-Weiss temperature scale. Neutron 
scattering reveals the presence of short range magnetic order- 
ing below 50K with an incommensurate wavevector; how- 
ever, the magnetic correlation length never grows past about 7 
lattice spacings even as the system is cooled to 1.5K.^ ~ More 
recent studies have shown that local moments do form and 
gradually freeze at low temperatures, fluctuating in creas ingly 
slowly as the temperature is reduced below They 
eventually become static on the scale of the slowest avail- 
able experimental local probe (muon spin resonance) below 



about 2K. 

A series of theoretical studies of this system subsequently 
appeared which attributed the thermodynamic signatures of 
ordering to a quadrupolar (also called spin-nematic) ground 
state stabilized by a n eares t-neighbor biquadratic coupling 
term in the Hamiltonian!^'^"^' Though such a state can account 
for the observed low temperature thermodynamics without 
requiring there to be any long-range spin ordering, it is in- 
consistent with the observed short-range incommensurate spi- 
ral order at the lowest temperatures. Moreover, the observa- 
tion of static (though spatially random) spins at low tempera- 
tures casts doubt on the relevance of the spin-nematic phase to 
NiGa2S4. Indeed, to date no clear specific signature of even 
quadrupolar correlations (much less order) has been identified 
in NiGa2S4. 

In this paper, we reconsider the principal experimental ob- 
servations, and argue that such correlations are likely present 
in NiGa2S4, and signified by the unusual two peak structure 
in the specific heat. More formally, the double peak is a signa- 
ture of proximity to a zero temperature quantum critical point 
between an antiferromagnetically ordered and a spin-nematic 
state. We also suggest a mechanism for the low temperature 
spin freezing, which is independent of the spin-nematic cor- 
relations. This mechanism should apply to any strongly two- 
dimensional magnet whose ground state in the absence of dis- 
order is an incommensurate spiral. We argue that both these 
phenomena (the two peak specific heat and spin freezing) can 
be understood on purely symmetry grounds. Thus they may 
be expected to occur much more generally in other materials 
under conditions that we discuss. 

For concreteness, we work within an extended model that 
not only includes the nearest-neighbor Heisenberg and bi- 
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quadratic couplings, but also third-nearest-neighbor Heisen- 
berg interactions: 



H 



M 



S,f. (1) 



The third neighbor couplings are denoted by [ij] (on the tri- 
angular lattice, the third neighbor bonds point along the first 
neighbor bonds but doubled in length). This additional J3 
term is perhaps the simplest way to account for the observed 
order, which is described by a wavevector that is slightly less 
than half of the 120-degree state wavevector. Indeed, such an 
incommensurate spiral state is the exact solution of the classi- 
cal limit of the above model in the limit J3 > —Ji > K > 0. 
Furthermore, detailed microscopic studies of NiGa2S4 indi- 
cate that a strong third-neighbor interaction is in fact present 
while second-neighbor interactions appear to be much less im- 
portant.^ 

Let us outline our approach to the problem, and the lay- 
out of the remainder of the paper. For J3 = 0, Eq. ([TJ has 
been analyzed already in Refs. I6I7I and shown to exhibit a 
non-magnetic but (ferro-)quadrupolar (spin-nematic) ground 
state in a broad parameter range. We argue that this remains 
true for the more general model with J3 ^ in Sec. |V] In 
addition, when K is not too large, Eq. ^ has magnetically 
ordered ground states. Each of these states is connected to 
the quadrupolar one by a quantum phase transition (QPT). A 
major portion of this paper is concerned with establishing the 
T > properties of systems in the vicinity of such a QPT. We 
first present general symmetry-based arguments for the exis- 
tence of two energy scales, and a resulting two peak structure 
of the specific heat, in Sec.|Il] To verify these ideas, we then 
introduce a novel "semiclassical SU(3)" or "sSU(3)" approx- 
imation, in which we can obtain these properties for Eq. ([T]i 
by classical Monte Carlo simulations. Though this approach is 
not expected to be quantitatively accurate for the 5 = 1 model 
of interest, it does preserve all the symmetries of Eq. (jTJ and 
a qualitatively appropriate T — limit, unlike conventional 
semiclassical approaches. We carried out extensive numerical 
simulations using the sSU(3) approach to establish the nature 
of the T > phase transitions and crossovers, which we de- 
scribe in Sees. HV] and [V] 

Next, we turn to the effects of disorder, assuming that the 
putative ground state in the absence of impurities is magnet- 
ically ordered (though it may be near the QPT separating it 
from the quadrupolar phase). We show that the generic in- 
commensurate spiral ground state of Eq. ([T]i is highly suscep- 
tible to defects. This is a direct consequence of the nature 
of the associated order parameter, which leads to strong "ran- 
dom field" effects of impurities - effects which do not pertain 
to the more conventional 120° three-sublattice ground state. 
Moreover, very general arguments due to Imry and Ma show 
that such random field systems in two dimensions - but not 
in three - are always disordered. Hence we are led to con- 
clude that the nature of the order parameter (which in turn 
is due to the large J3 interaction) and the very weak three- 
dimensional coupling in NiGa2S4 are responsible for the ob- 
served low temperature glassy state. This mechanism for spin 



a discussion in Sec. VII The Appendix describes details of the 
measurements used in the Monte Carlo simulations. 



II. SYMMETRY AND ORDER PARAMETERS 

In this section, we discuss the symmetries and order param- 
eters of the Hamiltonian in Eq. ([T]i and its various proposed 
ground states. The Hamiltonian itself has global SU(2) spin- 
rotation symmetry, and the full space group symmetry of the 
triangular lattice. It is also time reversal invariant. A subset 
of these symmetries are retained in the various ground states. 
First, we discuss the conventional magnetically ordered states. 
We consider only coplanar spirals, which dominate the classi- 
cal ground states of Eq. ([T]i. In these states, the pattern of spin 
ordering is described by 



(S,) = Rc [de'''-"-'] , 



(2) 



where d = \d\{ei + 162), with 61,62 orthogonal unit vec- 
tors, and k is the spiral wavevector. The spins lie in the plane 
spanned by ei, 62, normal to 63 = ei x 62. Clearly, the spin 
spiral breaks spin rotational symmetry, and, assuming k 7^ 0, 
it does so completely: rotation by any angle about any spin 
axis alters the spin state. Moreover, under the same assump- 
tion, the spiral breaks translational symmetry. However, an 
important symmetry is retained by the coplanar spiral. In par- 
ticular, under the combined action of a translation by any (lat- 
tice) vector a, and a spin rotation by the angle k • a about the 
63 spin axis, the state is unchanged. The situation for spatial 
rotations is more complex. Consider the three-fold (C3) rota- 
tion about a lattice site, a symmetry of the triangular lattice. In 
the familiar 120° state (the ground state for K = = 0), the 
ordered spins remain invariant under this rotation, since the 
sublattices are unchanged. This invariance is, however, not 
generic. Formally, it can be traced to the fact that in this case, 
k is a zone-corner wavevector (the "K point"), which is in- 
variant under such a rotation. Any other k is not invariant un- 
der a three-fold rotation, and moreover, the invariance cannot 
be regained by combining the spatial rotation with an SU(2) 
one. The breaking of the discrete C3 rotation symmetry of 
H implies that there are then three distinct "domains" of spi- 
ral state, with each of the three possible wavevectors forming 
an invariant set under the rotations. Each domain cannot be 
continuously transformed into another. This is quite different 
from the 120° state, in which all possible spin configurations 
can be smoothly transformed into one another. 

Next consider the (ferro)quadrupolar ground states. In them 
the ordered spin moment is zero: 



0. 



(3) 



However, the average magnetic quadrupole moment is non- 
zero: 



freezing is discussed in Sec. VI B We conclude the paper with 
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where is the (ferro-)quadrupolar order parameter By 
construction, Q is a traceless symmetric matrix. In the 
"spin-nematic" states which we consider, the eigenvalues 
of this matrix come in only two values, i.e. eigs{Q) = 
{q, ~q/2, ~q/2}. More complicated "biaxial spin-nematic" 
states are possible in principle, in which all three eigenvalues 
are distinct. We will not consider them here, as there is no 
obvious reason for them to occur in Eq. ([T]i. In the remainder 
of the paper, when we refer to the quadrupolar state, it will al- 
ways be assumed to be of the (single axis) spin-nematic type. 
Then a general Q can be written as 



(6) 



where n is a unit vector known as the "director". The change 
of the director by a sign, n — n, has no physical signifi- 
cance, because Q is even in n. The quadrupolar state breaks 
spin rotation symmetry, but retains invariance under spin rota- 
tions about the axis of the director. Furthermore, it retains full 
space group and time reversal symmetry. 

Comparison of the symmetries of the magnetically ordered 
and quadrupolar ground states reveals an important fact: the 
symmetry group of the spiral state is a subgroup of the sym- 
metry group of the quadrupolar one. Thus any order param- 
eter of the quadrupolar state must already be non-zero in the 
magnetic phase. Indeed, within Landau theory, one expects a 
coupling between the magnetic and quadrupolar order param- 
eters, i.e. a term in the free energy density / of the form 



•d,5, 



(7) 



where A is a non-zero coupling constant. This generically in- 
duces a non-zero Q^,^ of the form in Eq. (j6]), with \q\ cx jdp 
and n = §3. 

It follows that q in Eq. (j6]l remains non-zero across the QPT 
between the magnetic and quadrupolar states. Conversely, at 
this QPT, presuming it is continuous or weakly first order, the 
magnetic order parameter d becomes arbitrarily small. Asso- 
ciated with each of these order parameters is an energy scale, 
describing the energy needed to disturb the order parameter 
locally. This energy is supplied with increasing temperature, 
releasing entropy associated with the ordering. As usual, the 
specific heat is expected to show a peak at this characteristic 
temperature (whose precise nature requires more detailed con- 
siderations of symmetry, universality, and energetics). Near 
the QPT, on the magnetically ordered side, one is then led to 
expect two widely separated energy scales: the quadrupolar 
"ordering" scale, which remains of 0(1) (i.e. not parametri- 
cally small near the QPT), and a much smaller magnetic "or- 
dering" scale, which vanishes on approaching the QPT. There 
will be a specific heat peak at each of these temperatures. The 
schematic phase diagram is shown in Fig. [T] We have used 
"ordering" in quotes, because in a two-dimensional system, 
the continuous order parameters are prohibited from exhibit- 
ing long-range order by the Mermin- Wagner theorem. How- 
ever, the specific heat is dominated by short distance correla- 
tions, and exhibits peaks regardless of the presence or absence 
of true phase transitions. 
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FIG. 1: Finite temperature phase diagram of the nearest-neighbor 
only model with antiferromagnetic Ji > and a ferronematic K > 
obtained from numerical simulations in the sSU(3) approximation. 
The location of the phase boundaries are obtained from the position 
of the specific heat peaks in our numerical simulations. Approaching 
zero temperature the semiclassical sSU(3) approximation recovers 
the ground states obtained in mean-field theory. We note that the 
phase boundaries indicate a rapid growth of the ferroquadrupolar and 
magnetic correlation lengths, but not true long-range order excluded 
by the Mermin- Wagner theorem. 



III. SEMICLASSICAL SU(3) APPROXIMATION 

In the previous section, we argued for the emergence of 
two temperature scales and a two peak structure in the specific 
heat, whenever the system is in vicinity of the QPT between 
an antiferromagnetic and quadrupolar ground state. In what 
follows, we will check these arguments by explicit calcula- 
tions of T > properties. To do this presents a serious diffi- 
culty: theoretical approaches to calculate the low but nonzero 
temperature behavior of frustrated quantum Hamiltonians are 
extremely limited. A direct attack with quantum Monte Carlo 
is prevented by the sign problem, and many other approaches 
(e.g. variational wavefunctions, Lanczos diagonalization) are 
restricted to ground state properties. 

One natural approach is the classical approximation, in 
which the quantum spin operators are replaced by continuous 
classical vectors of length S. For many quantum spin models, 
the success of spin-wave theory, even for low spin, proves the 
applicability of this method. However, for the study of Eq. Q, 
such an approach is fatally flawed. The difficulty lies in the 
highly quantum nature of the quadrupolar spin-nematic corre- 
lations induced by the biquadratic interaction for spin S — 1. 
In particular, the ferronematic ground states favored hy K > 
are of "easy-plane" type, in which the spins may be thought 
of as fluctuating quantum mechanically in the plane perpen- 
dicular to the director h in Eq. (j6]l. This corresponds to the 
sign q < 0. For large spin, including the classical limit, how- 
ever, the biquadratic interaction with K > favors collinear 
states of the spins, and hence "easy axis" rather than easy- 
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plane ferronematic order (with the opposite sign of q). In a 
simple mean-field theory, the nature of the quadrupolar order 
changes from easy plane when 5* < 2 to easy axis for 5 > 2. 
Therefore even qualitatively the classical vector spin approxi- 
mation fails to capture the proper ferronematic correlations of 
the physics 5=1 system. 

We thus follow an alternative approach, which generalizes 
the zero temperature variational approach of Ref 7 to T > 0. 
In that work, a trial ground state was taken to be entanglement- 
free, i.e. of direct product form: 



I*) 



(g)i\'ipi 



(8) 



with an arbitary spin 5=1 state 1-0)^ on site i. A general 
single spin state can be written as 



= bii\x) + bi2\y) + bi3\z) , 



(9) 



with an arbitary complex vector satisfying the normaliza- 
tion constraint that b* • = 1. Here we have used the time- 
reversal invariant basis of 5 = 1 states defined by 



where |1),|0) and |1) are the usual 5=1 states quantized 
along the z-axis. These basis states are null vectors of their re- 
spective spin components, i.e. S^\x) — S'^\y) — S^\z) = 0. 
Semiclassically, the state |a;) may be viewed as one in which 
the spin fluctuates primarily in the yz-plane (and similarly for 
the other two states). 

The expectation values of the spin and quadrupolar opera- 
tors have simple expressions in these states: 



(V'dS.IV'^) = ~^b* xb, 



(10) 
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Note that when bj is real (it may be multiplied by an overall 
phase without affecting physical properties), the spin expec- 
tation vanishes. Thus in this case it may be interpreted as the 
nematic director, i.e. b^ = hi. Adding an imaginary part to 
hi produces a nonzero spin moment that is orthogonal to the 
plane in which the real and imaginary parts of b^ lie. 
From Eq. 



Hamiltonian, 



(11) 



J3^|b*.b,f -|b,.b, 



+ ^jiib:.b,f -(ji + if)ib,.b,f . 

m 

In Ref. |2l this expression was minimized for the nearest- 
neighbor model ( J3 = 0) to find variational ground states. 



It was found that the variational approximation qualitatively 
reproduced the correct T = phase diagram, though 
it incorrectly determined the quantitative location of the 
ferronematic-antiferromagnetic phase boundary. Note that, at 
T = 0, the variational approximation is equivalent to a mean- 
field theory in which the expectation values Qi-^f^, and (Si) are 
self-consistently determined from a Hamiltonian of decoupled 
5=1 spins. 

To address T > temperature properties, we follow an 
unusual procedure of elevating Eq. ([TTJ to a Hamiltonian of 
a two-dimensional statistical mechanical model, allowing the 
classical b^ vectors to fluctuate thermally with the Boltzmann 
weight exp[— /3iJci]- Because in this new classical problem, 
the "spins" appear as three-component complex vectors, and 
transform under symmetries by rotations by unitary three- 
dimensional matrices, we call this treatment the semiclassical 
SU(3) or sSU(3) approximation, and the resulting statistical 
mechanical problem the sSU(3) model. We emphasize, how- 
ever, that the model Eq. 



11 



does not have SU(3) symmetry, 
but only the physically appropriate SU(2) spin rotational and 
lattice symmetries. 

The sSU(3) model can be understood as a leading-order 
term in a cumulant expansion of the full quantum partition 
function. Specifically, we approximate 



Z = Tre 



-fJH 



(12) 



= J l[dn,,{^{b)\e-('"\^{b)) 



Y[da 



fc,e 



10 1, we can compute the expectation value of the Eq- P 1 



where drib = i'2^?[U^^^.,y,z ^R^^d^^WE^l 
1) is the appropriately normalized measure over unit 3- 
component complex vectors, such that J dfti,! = Trl = 3. 
From this perspective, we see that the sSU(3) approximation 
becomes exact at high temperatures (reproducing the exact 
free energy to 0{f3) in the p expansion), and approaches the 
variational result as T ^ 0. It can therefore be expected 
to yield a reasonable qualitative approximation to the physics 
over the full range of intermediate temperatures. 

One may consider at this stage, given that the T — vari- 
ational approach is equivalent to mean-field theory, why not 
instead apply the usual mean-field theory at T > 0? The rea- 
son to choose the sSU(3) approach is that it captures much 
better the effects of thermal fluctuations, which are very im- 
portant in a two-dimensional system. For instance, because 
has the physically correct symmetries, we expect by 



universality that it displays the correct exact critical behavior 
at any T > phase transitions. Moreover, it avoids spuri- 
ous transitions which are in fact prohibited by the Mermin- 
Wagner theorem. In contrast, straight mean-field theory will 
produce sharp phase transitions associated with quadrupolar 
and magnetic ordering, neither of which are in fact allowed in 
two dimensions. 

Despite the above advantages of the sSU(3) approximation 
over mean-field theory, it is very amenable to numerical com- 



putations. The classical Hamiltonian in Eq. ( 1 1 1 is completely 
real, so there is manifestly no sign problem. The only subtlety 
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lies in the gauge redundancy of the variables, whose phase 
(on each site) has no physical significance. However, no diffi- 
culty is incurred in Monte Carlo simulations by simply includ- 
ing the redundancy, and restricted measurements to gauge- 
invariant observables. For our numerical simulations of this 
semiclassical Hamiltonian we have built our implementation 
of a classical Monte Carlo algorithm on the ALPS libraries'^, 
and in particular expanded its classical Monte Carlo appli- 
cation code to sample configurations of the three-component 
complex vectors or "spins" as introduced above. We give a de- 
tailed account of how to measure various physical observables 
in this Monte Carlo scheme in the Appendix. 



IV. FINITE TEMPERATURE PROPERTIES OF THE 
NEAREST-NEIGHBOR MODEL 



As a first step we apply the semiclassical sSU(3) approx- 
imation to numerically analyze the T > thermodynamic 
features and correlations of the simplest, nearest-neighbor 
Hamiltonian with biquadratic exchange 



(a) 



H 



(13) 



where we consider antiferromagnetic Ji > and ferrone- 
matic K > 0. Though our primary focus will be on finite 
temperature properties of this model, it is helpful to begin by 
reviewing what is already known about its ground states. 



A. Ground State Properties 

The full zero temperature phase diagram of the nearest- 
neighbor bilinear-biquadratic model contains the conventional 
ferro- and antiferromagnetic magnetic phases when K is 
small, but also includes a ferroquadrupolar phase and an anti- 
ferroquadrupolar phase stablized by larger values of K. 

In the range of parameters of interest to us ( Ji and K pos- 
itive), there is a single quantum phase transition as shown in 
Fig.[T] For small K/ Ji, the ground state is a 120-degree spiral 
ground state with a reduced spin moment. Then, as K/Ji is 
increased, the spin moment decreases to zero leaving a state 
with only ferroquadrupolar order i.e. no average spin moment 
but an identical director (identical, real b) on every site. In 
the mean-field/variational approximation, this phase transition 
happens precisely at K/Ji = 2. Exact diagonalization re- 
sults^ for the quantum model show, however, that quantum 
fluctuations strongly reduce this value to about K/ Ji = 0.4. 

Within mean-field theory, a detailed study shows that the 
quadrupolar to antiferromagnetic quantum phase transitions 
is continuous. Thus this model fits nicely into the framework 
discussed in Sec. |ll] For any finite K/ Ji up to the critical 
value, the Neel and ferroquadrupolar orders actually coexist, 
leading to a smooth reduction in the value of | (S) | as K is 
increased. 
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FIG. 2: Finite temperature calculations of (a) the specific heat and 
(b) the spin and quadrupolar correlation lengths and for the 
nearest-neighbor model in the semiclassical sSU(3) approximation. 
The two temperature scales associated with the rapid growth of mag- 
netic and ferroquadrupolar correlations result in a characteristic two 
peak structure of the specific heat. The data shown is for couplings 
Kj ,]\ — 1.5 and system sizes chosen to be commensurate with the 
low-temperature spin ordering. 



B. Separation of Temperature Scales 

We have argued above that close to such a quantum phase 
transition from spkal to quadrupolar order a spin system 
may develop significant quadrupolar correlations at a distinct, 
higher temperature scale than the one associated with the rapid 
growth of magnetic correlations. Our numerical simulations 
of the sSU(3) model associated with Hamiltonian ( [T3] l does 
indeed reveal such a characteristic separation of temperature 
scales for K/ Ji sufficiently close to the quantum phase transi- 
tion which in mean-field approximation occurs at K/Ji — 2. 

In particular, we find that in the proximity of the quantum 
phase transition K/Ji < 2 the measured specific heat de- 
velops two distinct peaks as illustrated in Fig. |2] Calculat- 
ing the ferroquadrupolar and magnetic correlation lengths, we 
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can associate the upper peak with the rapid growths of fer- 
roquadrupolar correlations while the lower peak indicates the 
onset of magnetic correlations (with a wavevector of length 
47r/3 that is characteristic of the 120-degree state). We note 
that although the correlation lengths grow to be very large be- 
low each of the two temperature scales, they remain finite, 
consistent with the absence of long-range order While the 
high temperature peak in the specific heat shows no singular 
behavior as we cross the location of the quantum phase transi- 
tion at K/ Ji = 2, we find that the lower peak approaches 
zero-temperature exactly at the QPT. For smaller values of 
K/Ji < 2 the two peaks get closer and eventually merge 
around K/ Ji « 5/4, below which we only observe one tem- 
perature scale indicated by a single peak in the specific heat. 

We summarize our results for the nearest-neighbor model 
^T3\ in the phase diagram of Fig. [T] where the phase bound- 
aries indicate the locations of the respective peaks in the cal- 
culated specific heat. Note that these phase boundaries are not 
phase transitions, due to the absence of long-range order as 
mandated by the Mermin- Wagner theorem. 

Finally, we note that the numerical results of our sSU(3) 
approximation are not expected to be quantitatively correct, 
as it becomes clear from the quantitative discrepancy to ex- 
act diagonalization results. However, we emphasize that the 
key ingredient underlying our qualitative results is a contin- 
uous quantum phase transition in which the magnetic order 
smoothly vanishes but the quadrupolar order survives in both 
phases. As long as the full inclusion of quantum effects does 
not lead to the zero temperature phase transition actually be- 
coming strongly first order, then a separation of temperature 
scales should generically be observed when the parameters of 
this system are near the quantum critical region. 
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V. THE J3-J1-K MODEL 

We now turn to the full model in Eq. ([T]) including a third- 
nearest-neighbor coupling, which is necessary to properly 
capture incommensurate magnetic correlations at low temper- 
atures. We assume antiferromagnetic third-neighbor interac- 
tions J3 > and ferronematic K > {as before), but now 
consider a /erromagnetic nearest-neighbor coupling Ji < 0. 
It is this combination of couplings that will realize the proper 
magnetic correlations observed in NiGa2S4. As we shall see, 
the presence of an incommensurate spiral ground state has in- 
teresting consequences in that it introduces a true finite tem- 
perature phase transition associated with anisotropic spin fluc- 
tuations. 

While we find that the above model can now account for 
the type of magnetic order observed in NiGa2S4, it still can- 
not explain the fact that the spin order never grows to be long 
ranged and exhibits slow dynamics. Therefore we will ulti- 
mately have to consider the role of disorder before discussing 
our results in the context of the real material. 
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FIG. 3: Zero temperature phase diagram of the model with antifer- 
romagnetic third-nearest-neighbor coupling J3 and nearest-neighbor 
ferromagnetic Ji and ferronematic K couplings obtained from nu- 
merical simulations in the sSU(3) approximation. The three phases 
present at zero temperature have a ferroquadrupolar, 120 degree spi- 
ral or incommensurate spiral ground state, with continuous transi- 
tions (dashed lines) between the ferroquadrupolar and spiral states 
and a first-order transition between the 120 degree and incommensu- 
rate spiral states. The cross-hatched area marks the region where our 
finite temperature calculations indicate two temperature scales. Such 
a separation is observed both in specific heat (filled diamonds) and 
correlation length measurements (open circles) which are expected 
to coincide in the thermodynamic limit. The red triangle marks the 
location of the representative system, for which we show finite tem- 
perature properties in Fig.|4] 



A. Ground States 

We again start our analysis of of Hamiltonian ([T]) by dis- 
cussing its zero temperature phase diagram calculated in 
mean-field approximation. For the range of parameters of 
interest, we assume that it is sufficient to take the magnetic 
ground states to be coplanar and invariant under the combina- 
tion of a lattice translation and spin rotation associated with 
a single wavevector. Under these two assumptions, the most 
general possible choice for the vector of amplitudes de- 
scribing the wavefunction at a single site is, up to a global 
0(3) spin rotation. 



sin 6* cos (k • ri)e' 
sin6'sin(k • ri)e" 
cos 6 



Minimizing the energy over the variables k, 6 and a we obtain 
the zero temperature phase diagram shown in Fig. [3] For small 
values of K/ J3, the ground states are again spirals with a re- 
duced spin moment. Along the line | Ji| = K/2 for K < 1, 
the wavenumber describing these spiral states ground states 
jumps from 47r/3 to 27r/3 and then continuously decreases to 
zero with increasing |Ji|, ultimately leading to a feiTomag- 
netic ground state for very large values of | J1I/J3. The new 
feature in this phase diagram is thus a large region of parame- 
ters for which the spiral is incommensurate (with wavevector 
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k < 27r/3). Increasing the biquadratic exchange K for fixed 
\Ji\/J3^ these spiral ground states smoothly transform into 
a ferroquadrupolar ground state with (S) = which is the 
ground state for large K. 



B. Separation of Temperature Scales 

Similar to the nearest-neighbor only model we expect on 
symmetry grounds that also the model with third-neighbor in- 
teractions exhibits a separation of temperature scales in the 
vicinity of the continuous quantum phase transition between 
the ferroquadrupolar and (incommensurate) spin spiral state. 
Indeed our numerical simulations of this model in the sSU(3) 
approximation reveal a broad range of couplings in the vicin- 
ity of this line of continuous transtions where both specific 
heat and correlation length measurements indicate two differ- 
ent temperature scales as indicated by the cross-hatched area 
inFig.||] 

In agreement with the symmetry arguments in support of 
this temperature scale separation we again find that the up- 
per temperature scale corresponds to the rapid growth of 
quadrupolar correlations while the lower temperature is asso- 
ciated with the onset of magnetic correlations. Numerical data 
for a representative choice of parameters (indicated by the red 
triangle in Fig. [3]) is shown in Fig. |4] where we plot the spe- 
cific heat and correlation length measurements as before. The 
exact parameters K/J^ — 1.1 and J1/J3 ~ —1.14855 and 
system sizes are again chosen such that the low temperature 
spiral state has a commensurate spatial period (in this case 7 
lattice spacings) which in turn allows to compute the correla- 
tion lengths directly from the respective structure factors. 



C. C3 Bond Ordering Transition 

As ferroquadrupolar correlations only break the continuous 
SU(2) symmetry, the transition at the upper temperature scale 
should not be associated with the onset of true long-range or- 
der as mandated by the Mermin-Wagner theorem. Indeed a 
finite-size scaling analysis of our numerical data shows that 
the high-temperature peak in the specific heat is non-divergent 
as we increase the system size, as illustrated in Fig.|5] 

A similar analysis for the transition at the lower temperature 
scale, however, reveals that more interesting things can hap- 
pen. For any choice of parameters that puts the ground state 
of the system into the incommensurate spiral phase, the lower 
peak of the specific heat appears to diverge sharply as shown 
in Figs. |4|a) and [6] This divergence indicates a true phase 
transition, which in two dimensions must reflect the breaking 
of some discrete symmetry in accordance with the Mermin- 
Wagner theorem. (In principle such a divergence could also 
be explained by a first-order transition without a change of 
symmetry, but will be ruled out below.) 

The key observation to explain this phase transition to a 
long range ordered state is that any (incommensurate) spiral 
state on the triangular lattice with wavenumber |fc| < An/ 3 
exhibits spin correlations that are stronger along one axis than 
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FIG. 4: Finite temperature calculations of (a) the specific heat, 
(b) the spin and quadrupolar coiTelation lengths and and (c) 
the nematic order parameter \ip\'^ for the model with third-nearest- 
neghbor interaction in the semiclassical sSU(3) approximation. The 
two temperature scales associated with the rapid growth of mag- 
netic and ferroquadrupolar correlations again result in a character- 
istic two peak structure of the specific heat. The data shown is for 
couplings K/Js — 1.1 and J1/J3 = —1.14855 and system sizes 
chosen to be commensurate with the low-temperature spin ordering. 
In contrast to the nearest-neighbor model, we see a divergent lower 
peak corresponding to the breaking of lattice rotational symmetry by 
anisotropic spin fluctuations. 
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FIG. 5: Finite-size scaling of the upper peak in the specific heat of 
the representative system with couplings Kj Jz = 1.1 and ,h/ Jz = 
— 1.14855. While the peak slightly shifts with increasing system 
size, there is clearly no divergence. 




0.225 0.23 0.235 0.24 0.245 0.25 0.255 



temperature T / 

FIG. 6: Finite-size scaling of the lower peak in the specific heat 
of the representative system with KjJz ~ 1.1 and J1/J3 ~ 
— 1.14855. The trend with increasing system size clearly suggests 
a divergence in the thermodynamic limit. 



the other two, unlike the 120-degree state, in which every spin 
differs from its neighbors by the same angle. Therefore any 
incommensurate spiral state spontaneously breaks the discrete 
C3 rotational symmetry of the lattice and there is nothing to 
prevent this symmetry breaking from occurring at a finite tem- 
perature. Indeed, it has recently been shown that the related 
classical model with 0(3) spins on the triangular lattice and 
first and third neighbor exchange spontaneously breaks the C3 
lattice rotation symmetry at low temperature'^. In the C3 bro- 
ken state the system exhibits "bond order" in the sense that 
the exchange energies for bonds along different principle axes 
become distinct. 

To confirm that the lower temperature scale in our sSU(3) 
model does indeed correspond to the breaking of lattice rota- 
tional symmetry, we define a bond order parameter ip, which is 




-2.890 -2.878 -2.866 -2.854 -2.842 -2.830 

2 

energy density E / (J^ L ) 

FIG. 7: Histogram of the energy density for our representative 
system in the vicinity of the low temperature transition. The bi- 
modal character of this distribution proliferates with increasing sys- 
tem size, which is indicative of a first-order transition. Data for dif- 
ferent system sizes/temperatures are shown; in particular L = 84 
(T/Ji = 0.2377), L = 98 (T/J3 = 0.2374), and L = 112 
(T/J3 = 0.2361). 



nonzero only in the presence of anisotropic spin correlations: 

where ai — x and 3.2.3 are its rotations by 2tt/3 and 47r/3, re- 
spectively. This order parameter is a complex number whose 
argument reveals the direction of the strongest correlations 
and whose magnitude for a spiral state with wavenumber k 
is given by 

IV'I = (|S|)2 [cos(fc) - cos(A:/2)] . 

Moreover, we note that ip is invariant under the transformation 
Si —Si and under a lattice rotation of angle tt. Thus, any 
nonzero value of tp indicates that the system is one of three 
inequivalent bond ordered phases related to each other by lat- 
tice rotations of 2tt/3 (cf. the discussion on spiral magnetic 
domains in section|Il]). 

Measurements of this bond order parameter for the sSU(3) 
model associated with Hamiltonian ([T]) confirm that the sys- 
tem enters a bond ordered state for couplings leading to an 
incommensurate spiral ground state, as shown in Fig.|4](c) for 
our representative system. To determine the order of this low 
temperature transition we have measured energy histograms 
in the vicinity of the transition temperature. The bimodal his- 
togram distribution, shown in Fig.|7j unambiguously points to 
a first-order transition, which is also indicated by the naiTow 
peak in the specific heat and the apparent jump of the mag- 
netic coiTelation length in Fig.[4](b). Thus the C3 transition 
in the sSU(3) model has the same qualitative properties as the 
corresponding transition in the 0(3) model which was also 
found to be first order'^. 
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VI. DISORDER 

A. Imry-Ma argument and correlation lengths 

We now turn to the effects of disorder. We consider sul- 
fur (S) vacancies to likely be the dominant type of impurity 
in NiGa2S4. However, our analysis will be rather general, 
and rely only upon the assumption that the defects are un- 
correlated, local, and non-magnetic. On symmetry grounds, 
non-magnetic disorder couples directly to the C3 order pa- 
rameter, but not to d, since (neglecting spin-orbit) it respects 
SU(2) symmetry. Thus for the C3 order, the impurities act 
as "random fields". In two dimensions, the standard Imry- 
Ma argument implies that long-range C3 order is destroyed 
by arbitrarily weak disorder.ElThe system breaks up into "do- 
mains", leading to a saturation of the correlation length 
for the C3 bond order. We parametrize the strength of the 
random field by an energy A, which gives e.g. the range of 
variation in strength of exchange couplings due to S vacan- 
cies. The essence of the Imry-Ma argument is a competition 
between the surface energy required to nucleate a domain of 
one phase within another, and the energy gain that can be ob- 
tained from the random fields within the domain. The former 
is proportional to the surface tension cr of a domain, while the 
latter is proportional to A 1-01. Hence we expect to be a 
function of a/(A|'0|). According to the Imry-Ma argument,^^ 
two dimensions is the marginal (lower critical) dimension for 
the stability of the ordered state. Hence the correlation length 
is an exponential function of the disorder strength. We arrive 
at the conclusion 



A2|0(T)|2 



(14) 



with some constant c. At T ^ Tc,. we expect a{0) ~ 0(^3) 
and ■!/'(0) ^ 0(1). Eq. (14i describes a continuous decrease 



of the C3 correlation length on increasing temperature up to 
T Ri Tc^, the temperature of the C3 transition of the pure 
system. At this point, the first order transition of the pure 
system will be rounded by disorder, and the C3 correlation 
length will smoothly (and quickly) approach its finite value 
for the pure system for T > Tc.^ ■ 

What of the magnetism? Although the magnetic order 
parameter does not couple direction to disorder, it will be 
strongly affected by the domain structure. In particular, the 
conventional spkal correlation length will saturate at the do- 
main size. 



6/(A,T) «min [^m{Q,T),£,cA^,T)] . 



(15) 



Because the pure system's magnetic correlation length grows 
very large below Tc^ , we expect that for most if not all of this 
temperature range, that simply = = ^■ 



B. Correlations on longer scales: gauge glass physics 

The finiteness of the above correlation lengths does not im- 
ply a complete absence of magnetic correlations beyond this 



length scale, only that the spiral correlations remain short 
range. One longer scales, we need to consider more subtle 
correlations between spins in different domains. 

In what follows, we will make a simplifying assumption, 
that the spirals are fixed in the XY plane. This could be justi- 
fied by a small easy-plane anisotropy, i/anis = ^X^i {^f)^^ 
which is expected to be present, but small, in NiGa2S4. Pro- 
vided Z?^^ > /csT, which will be satisfied at low enough 
temperatures, the spins will indeed lie in the XY plane. It may 
also be the case that the spins adopt a coplanar state sponta- 
neously, even in the absence of magnetic anisotropy. How- 
ever, we do not make a definitive statement in this respect. 
Should this assumption fail, we expect that this approximation 
underestimates the fluctuations, and therefore overestimates 
the tendency to spin freezing. One should probably therefore 
view the estimates of the spin glass correlation length and time 
below as upper bounds. 

To proceed, we define = (x + iy)e'^' as the average or- 
der parameter for a domain i, within which the wavevector is 
fixed at Qb.. In general, the energy of the system will depend 
upon the difference of the angles 6i — 9j in adjacent domains, 
due to the coupling of spins at the domain boundaries. This 
coupling is complicated and depends upon the specific struc- 
ture of this boundary, and the wavevectors of the two domains. 
Therefore we will assume the effective Hamiltonian 



eff 



cos 



(16) 



where aij is a random variable, which we take to be indepen- 
dent on each ij and uniformly distributed between and 27r. 
We have neglected for simplicity randomness in the magni- 
tude of the coupling between domains, which does not qual- 
itatively affect the physics. The coefficient K arises from 
adding the exchange couplings along the boundary of the two 
domains in question. If the addition is constructive, we ex- 
pect, ignoring fluctuations K ~ Js^, or if it is random in- 
stead K ^ JsVC- With thermal and/or quantum fluctua- 
tions, this should be reduced by a factor of Assum- 
ing a mean-field temperature dependence of the magnitude of 
the local spiral order parameter, one finds a smooth reduction 
with increasing temperature K{T) = if (0)(1 - T/Tmf), for 
T < Tmp. Generally, K ^ fc^T in the regime where we can 



apply Eq. ([16 1. 

Eq. ( [T6] l is the Hamiltonian of the two-dimensional gauge 
glass model, which has been intensively studied in the 
past-'StSmni It is believed that its lower critical dimension 
is greater than two, so that the system is disordered (param- 
agnetic) at all T > However, the correlation length for 
glassy (Edwards-Anderson) order, and the correlation times 
of the spins, diverge rapidly as T ^ 0. 

To understand this behavior in more detail, we must dis- 
cuss the excitations of the gauge glass ground state. Obvi- 
ously, because of the global U{1) symmetry of Eq. ( 16 1, there 



are Goldstone-type spin waves, in which the phases 6i vary 
slowly away from their ground state values. These are the 
Halperin-Saslow modes, introduced in Ref . [Tsl and discussed 
in Ref. 19 in the context of NiGa2S4. More subtle, however, 
are non-smooth deformations of the phase, which have lower 
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energy than the spin waves for a corresponding length scale. 
These are the "droplets" of the droplet theory. From numer- 
ical studies, these have a classical energy E - AL-^^K with 
6 Ki —0.36.'^ Here L should be measured in units of the 



effective lattice spacing of Eq. (16 1, and A should be of order 
K. Thus at low temperature, those "droplets" with size larger 
than the glass correlation length. 



(17) 



will be thermally activated, and presumably fluctuating in and 
out. Here fc^To ^ K. These droplets should be responsible 
for the vanishing of the static moment. 

The question is, what is the time scale for these fluctu- 
ations? This is generally determined by the energy barrier 
which must be overcome to create one of the droplet excita- 
tions (these can also be thought of in terms of moving vortices 
in some region of size ^g). According to scaling theory, this is 
of order B£^^, with some "barrier exponent" tp. This gives the 
general expectation 



exp 



(To/T) 



i+i,/\e\ 



(18) 



Some recent numeric^ suggest t/j — 0, so that the barriers 
are logarithmic, of order B In (again B ^ K). In this case, 
the glass correlation time behaves like 



Tg ~ explBln^g/ksT] exp 



Toln(ro/T) 
T 



(19) 



Both Eq. ([T8| and Eq. ( [T9] l describe super-Arrhenius diver- 
gence of the relaxation time at low temperatures. For times 
shorter than this correlation time, the spins are effectively 
frozen. Due to the rapid growth of the correlation time at low 
temperature, this gradual process can be easily misinterpreted 
as a sharp dynamical transition. Nevertheless the actual freez- 
ing within the gauge glass model is a gradual and smooth pro- 
cess. Specifically, dependent upon the probe used, we expect 
that persistent spin dynamics can be observed below a nomi- 
nal "freezing" temperature (most likely the latter, empirically 
determined, temperature will be close to Tq^). 



sensitive to disorder, which acts like a "random field" on the 
C3 order parameter We argued that this results in a gradual 
freezing phenomena, related to the well-studied "gauge glass" 
model, with persistent but rapidly slowing dynamics down to 
zero temperature. 



B. Relation to experiments 

1. Model 

Having established these theoretical results, we now dis- 
cuss their relevance to NiGa2S4. First, we consider the ap- 
propriateness of the model Hamiltonian, Eq. ([T]l. Neutron 
experiments have decisively established the nature of short- 
range spin correlations. Their location near the ^,0) point 
in the Brillouin zone clearly implies the dominance of third 
neighbor antiferromagnetic exchange J3 over first and second 
neighbor exchange. Classically, weak ferromagnetic Ji se- 
lects this ordering wavevector and accounts for the slightly 
incommensurate value. This is also true quantum mechani- 
cally: for a model involving only Ji and J3, with Ji ^ J3, it 
is possible to confirm this type of spkal ground state directly 
for S = 1, assuming only that a nearest-neighbor triangu- 
lar antiferromagnet with S' = 1 is ordered.*^ This pattern of 
exchange interactions has also been justified by microscopic 
considerations.'^'^ On physical grounds, this situation implies 
an anomalously low overlap in the nearest-neighbor Ni-S-Ni 
superexchange pathway, which involves a 90° bond. Accord- 
ing to the usual Goodenough-Kanamori rules, such a pathway 
indeed results in a weak ferromagnetic exchange. However, 
this situation is very sensitive to variations in the bond an- 
gle. We may therefore expect virtual fluctuations of phonons 
which deform this bond angle to make an important contribu- 
tion to the effective magnetic Hamiltonian. In the limit of a 
high energy phonon (with /ic^phonon ?J Ji), this directly leads 
to the effective ferro-biquadratic coupling K in Eq. ([T}. 



2. Intermediate energy physics 



VII. DISCUSSION 
A. Summary 

In this paper, we have studied an extended bilinear- 
biquadratic exchange model for an 5 = 1 antiferromagnet on 
the triangular lattice, focusing on T > properties. We found 
that, when the system is proximate to a T = quantum phase 
transition to a ferro-quadrupolar (spin-nematic) state, the sys- 
tem is characterized by two energy scales. The specific heat 
then displays separate peaks at the two corresponding tem- 
peratures. We introduced a novel sSU(3) approximation and 
associated numerical method to study such T > proper- 
ties. Below the lower temperature peak, long-range C3 bond 
order is present, and spiral magnetic correlations rapidly de- 
velop. In this range of temperature, the system is exquisitely 



Is there any empirical evidence for interactions, such as K, 
beyond Heisenberg exchange? Though much of the empha- 
sis, both theoretically and experimentally, on NiGa2S4 has 
been on its low temperature (below say lOK) behavior, the 
high temperature phenomena (say for lOK < T < lOOK) is 
perhaps more anomalous. In particular, the two peak struc- 
ture of the magnetic specific heat, with no indications of a 
sharp phase transition at either temperature, is quite unusual. 
Moreover, the very high temperature of the upper peak, com- 
parable to the Curie- Weiss temperature itself, is especially 
striking for a frustrated magnet, in which the entropy usually 
remains down to temperatures well below Qcw- Moreover, 
Ga nuclear magnetic resonance (NMR) and muon spin res- 
onance (/iSR) experiments also show strong temperature de- 
pendence of the dynamics in this broad temperature range. In 
our opinion, these observations strongly suggest that a simple 
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picture in which only "classical" frustration is operative is in- 
adequate. Indeed, it is natural, given the coincidence of the 
Curie-Weiss temperature and the upper specific heat peak, to 
think that some unfrustrated type of correlations are at work. 
The ferro-quadrupolar order parameter discussed here is our 
proposed candidate, and we take the two-peak structure of the 
specific heat as an indication of the proximity of the ground 
state to a T quantum phase transition to a non-magnetic 
quadrupolar phase. This is fully consistent with the presence 
of non-negligible K > 0. 

It would be highly desirable to elaborate more quantita- 
tively upon this scenario. It is tempting to apply directly the 
Monte Carlo results of this paper within the sSU(3) approx- 
imation to fit data on NiGa2S4. However, we believe this is 
unjustified, because the sSU(3) approximation is known to be 
quantitatively bad. As shown in Ref. i at T = it dras- 
tically overestimates the biquadratic exchange K required to 
stabilize the ferro-quadrupolar phase in the nearest-neighbor 
model. Significantly more theoretical work is required to im- 
prove on the sSU(3) approximation for the extended exchange 
model appropriate for NiGa2S4. The low frequency dissipa- 
tive susceptibility, relevant to the NMR and /xSR 

experiments, is particularly worthy of study. 



3. Low temperature freezing 

Experiments have uncovered a complex freezing process 
in NiGa2S4. The temperature Tf « IQK is identified as a 
"freezing temperature" at which there is an apparent diver- 
gence of the l/Tjand 1 /T2 relaxation rates in NMR and NQR 
measure me nts.EEEl However, static moments are not formed at 
this temperature, and several measurements indicate that spin 
dynamics persists down to temperatures of at least 2K. For 
T < 2K, broad NQR and NMR spectra are observed, as is 
a 1/3 "tail" in the /iSR spectrum; both are indicators of in- 
homogeneous moments that are static on times of at least the 
(slowest) muon lifetime of w 2.2/isec. Relaxation of the 
1/3 tail is observed for T > 2K, indicating that here the spin 
dynamics is faster than r^. This is consistent with the lack of 
an NQR signal for 2K < T < Tj, which indicates that the 
spin dynamics has substantial weight at the NQR frequency. 

In our model, the low temperature specific heat peak is as- 
sociated with the onset of substantial magnetic correlations 
and C3 bond order. Experimentally, this coincides approxi- 
mately with Tf w lOK. Around this temperature, we ex- 
pect an increase of the magnetic correlations towards their low 
temperature saturation value determined by the Imry-Ma dis- 
order physics. As temperature is lowered below this value, we 
cross into the gauge glass regime. Therefore one expects in 
this range a gradual freezing process with spin dynamics gov- 
erned by Eqs. ( 18|19 1. This growth of glassy relaxation time 



Tg appears consistent with the experimental observations. In 
this picture, there is no true transition at the experimental Tf. 
Rather, what is observed experimentally are the residual ef- 
fects of the true transition that would occur around this tem- 
perature in an ideal clean sample. 



4. WhyNiGa2Si? 

A number of aspects of the experimental behavior of 
NiGa2S4 appear rather unique among frustrated magnets, in- 
cluding a variety of other triangular lattice systems which have 
been studied. The presence of the double peak specific heat, 
the absence of any magnetic or structural transitions, and the 
gradual nature of the freezing, are all unusual. It is therefore 
interesting to ask what makes NiGa2S4 special? 

Very likely there are several features of the material which 
conspire to produce this physics. First, the unusual exchange 
paths, which are a result of the Ni-S geometry, are critical 
for weakening the nearest-neighbor exchange. This allows for 
relatively large third-neighbor exchange and biquadratic cou- 
pling. The latter can probably understood as arising because 
although the microscopic overlaps of Ni d-orbitals with the 
S p-states are not small, there is an accident of bond angles 
which is responsible for the small Ji exchange. As a result, 
fluctuations of this angle lead to a large response, and inte- 
grating the corresponding phonons out naturally leads to sub- 
stantial nearest-neighbor biquadratic coupling of the sign pos- 
tulated here. The large J3 exchange leads to the incommen- 



surate correlations, which in turn we have argued in Sec. VI A 



are responsible for the lack of magnetic order and spin freez- 
ing. The large K is responsible for the quadrupolar correla- 
tions and hence the double peak specific heat. 

Also crucial to the analysis was the two dimensionality of 
the system. In general, there is also some interlayer coupling. 
But the very large separation between Ni layers in NiGa2S4 
makes this exceedingly weak. The nearly complete absence 
of correlations between layers has indeed been verified in neu- 
tron scattering.^ Substantial interlayer coupling could lead to 
sharp magnetic and/or quadrupolar transitions, which are not 
observed. 

Finally, a high degree of magnetic isotropy is required to 
justify the Heisenberg treatment of the system. Sufficient 
magnetic anisotropy would render the material XY-like, lead- 
ing to Kosterlitz-Thouless type physics, and also, if strong 
enough, further smoothing the quadrupolar ordering crossover 
(since explicit magnetic anisotropy acts like a symmetry- 
breaking field to the quadrupolar order parameter). Due to 
the closed shell configuration of the Ni^+ ions, the single- 
ion anisotropy is indeed expected to be tiny. Indeed, very 
small anisot ropi es have been inferred from ESR and NQR 
experiments.'^ 



[3l2n 



C. Relation to prior work 



Prior theoretical worlP^H pointed out the possibility of 
spin-nematic phases on the triangular lattice, and suggested 
these as candidate ground states for NiGa2S4. Though this is 
clearly not the case, the present paper builds upon these works 
and preserves quadrupolar correlations as an essential ingre- 
dient. 

More recently, there have been two Monte Carlo investiga- 
tions of classical 0(3) spin models related to ours. Tamura and 
Kawashima studied the model with first and third neighbor 



12 



Heisenberg exchange, but no biquadratic coupling." 'They ob- 
served a first order transition into a low temperature bond or- 
dered phase with broken C3 symmetry, similar to what occurs 
in our model. Kawamura and Yamamoto studied the model 
with nearest-neighbor Heisenberg and biquadratic exchangdSl 
(similar to what we discuss in Sec. |IV[ but crucially different 
in that they use 0(3) spins rather than the sSU(3) approach). 
They indeed observe two temperature scales, with the upper 
corresponding to quadrupolar correlations and the lower to 
spin correlations. However, because of the classical spin treat- 
ment, the nature of these correlations is qualitatively different 
from what occurs for spin 5 = 1 (and what we observe here). 
In particular, in their analysis the quadrupolar correlations are 
of easy axis type rather than easy plane, and the magnetic cor- 
relations are thus coUinear rather than spiral. 



broadly to many magnetic materials. It should be interesting 
to seek applications of these ideas amongst other frustrated 
magnetic materials. 
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APPENDIX: DEFINITION AND IMPLEMENTATION OF 
MEASUREMENTS 



D. Outlook 



1. Theoretical issues 



This work raises a number of theoretical problems which 
might be addressed in the future. Our sSU(3) treatment cap- 
tures the essential quantum nature of the phases, but does not 
properly treat the quantum critical point (QCP) itself The 
field theory of the QCP is interesting, insofar as it may be im- 
portant to consider the coupling of the gapless "quadrupole 
waves" (Goldstone modes of the spin-nematic order) to the 
magnetic order parameter In turn, we may expect the uni- 
versal properties of this QCP to govern the physical behavior 
in the temperature range between the two specific heat peaks. 
Indeed, even classically, the spin dynamics in this temperature 
regime was not addressed here, and is important to understand 
for comparison to a considerable amount of NQR data. 

It would be highly desirable to improve upon the sSU(3) 
approximation, to bring it into better quantitative agreement 
with exact diagonalization results at T = 0. At T = 0, one 
possible approach would be to adopt a description in terms of 
"tensor product states". The direct product variational form 
to which the sSU(3) approximation reduces at T = can 
then be viewed as the simplest, zero-entanglement type of 
such a tensor product state. Increasing the size of the tensor, 
and thereby introducing some entanglement, would systemat- 
ically improve the ground state wavefunctions. Generalizing 
the sSU(3) approximation to extend the tensor product state 
description beyond zero temperature is an interesting theoret- 
ical challenge. 



2. Generality of the results 

Though much of the paper was devoted to a Monte Carlo 
study of the Ji~ J^ — K model on the triangular lattice within 
the sSU(3) approximation, the arguments presented in Sees. [Ill 
and VIA are quite general. In particular, the emergence of 
two temperature scales near a quadrupolar to magnetic QPT, 
and the sensitivity of two-dimensional incommensurate spi- 
ral states to non-magnetic impurities should both apply rather 



In the following we give a detailed account of how to re- 
late measurements of physical observables in the original spin 
models ([TJ and ([13]) to expressions in terms of the three- 
component complex vectors of our semiclassical sSU(3) ap- 
proximation. 

We implemented measurements of the spin structure factor 
iSk by making use of the identification 

Si = -ih* X h, . 

We then computed the structure factor by taking 

A. = E< 



5k = (|Akp) . 

Similarly, we computed the quadrupolar structure factor Qk 
by measuring the quadrupolar order parameter defined by 



1 



1 



and then computing 



2k = Ed^ri 



Having obtained the spin and quadrupolar structure fac- 
tors, the spin and quadrupolar correlation lengths for ordering 
with wavenumber k were obtained by assuming the Ornstein- 
Zernike (mean-field) form for the structure factor peaks. One 
then obtains the expressions 




Here (Sp) stands for an average of the structure factor over 
wavevectors of magnitude p oriented along each of the six tri- 
angular lattice bond directions. The shift S can be any small 



13 



number for the case of the infinite system, but for the triangu- 
lar lattice of linear dimension L, we took 6 — 'iir/L since the 
structure factors are well defined only at the discrete Born- 
von Karman momenta. Finally, since we assume above that 



the structure factors are mean-field like, the results for the 
correlation lengths can only be trusted above the respective 
temperature scales at which each correlation length grows to 
be large. 
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